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' The neutron pair correlation in nuclei near the neutron drip-line is investigated using the selfcon- 
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sistent continuum Skyrme-Hartree-Fock-Bogoliubov theory formulated with the coordinate-space 
Green's function technique. Numerical analysis is performed for even-even = 86 isotones in 
the Mo-Sn region, where the 3p^/2 ^-^id 3pi/2 orbits lying near the Fermi energy are either weakly 
bound or unbound. The quasiparticle states originating from the I = 1 orbits form resonances with 
large widths, which are due to the low barrier height and the strong continuum coupling caused 
by the pair potential. Analyzing in detail the pairing properties and roles of the quasiparticle 
resonances, we found that the / = 1 broad quasiparticle resonances persist to feel the pair potential 
and contribute to the pair correlation even when their widths are comparable with the resonance 
energy. 
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I. INTRODUCTION 



Many interests have been taken recently on pairing properties of neutron-rich nuclei near 
the drip-line. The most peculiar case could be the firstly observed halo nucleus, ^^Li, where 
two neutrons forming the halo would not be bound to the nucleus if the pair correlation 
were absent [l-4|. Similar examples are also suggested in nuclei near the neutron drip- 
line in heavier mass region: e.g., possible giant halo structure (consisting of more than 
two neutrons) which is predicted in Ca and Zr isotopes by the selfconsistent mean-field 
models 5|-|7|. A new aspect in these examples is that the pair correlation occurs among 
neutrons occupying unbound or weakly bound orbits with low angular momentum / = 
or 1 whose wave functions extend far outside the nucleus due to the low (zero) centrifugal 
barrier. 

A useful theoretical tool to study the pair correlation in the weakly bound nuclei in all 
the mass regions, except the lightest ones, is the coordinate-space Hartree-Fock-Bogoliubov 
(HFB) approach js-lll, in which the quasiparticle wave functions of weakly bound and 
unbound nucleons are described in the coordinate-space representation. Indeed the pair- 
ing properties in nuclei near the drip-line have been studied extensively within the HFB 



scneme 



HQ 



21| as well as the relativistic Hartree-(Fock)-Bogoliubov models j4-|6| 



22 



25| . The pair correlation we have to deal with here is, however, a rather complex and un- 



resolved problem, and there exist controversial issues concerning the roles of weakly bound 



and unbound orbits. For instance, it has been argued in Refs. |18l-l2l| that neutrons in the 
weakly bound and unbound orbits with / = and 1 tend to decouple from the pair field gen- 
erated by the other neutrons because of the large spatial extension of their wave functions. 
It is also claimed that those neutrons contribute very little to the total pair correlation in 
nuclei. On the contrary, other studies show large pairing effects even on the weakly bound 
neutrons, leading to the pairing anti-halo effect 26| and the increase of the neutron pairing 
gap for weaker binding of neutrons or for shallower neutron Fermi energy 



In this paper, we would like to present an investigation of the pairing properties of nuclei 
close to the neutron drip-line, with intentions to clarify the roles of weakly bound and 
unbound orbits with low angular momenta. 

To perform this study, there exist some physically and technically important points which 
need to be treated carefully. Firstly, precise description of the wave functions outside the 
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nucleus must be guaranteed since we deal with weakly-bound and unbound orbits. We 
achieve it in the present study by using the coordinate-space mesh representation for the 
Skyrme-Hartree-Fock-Bogoliubov model 

Secondly, also related to the first point, a suitable boundary condition needs to be im- 
posed on the wave functions of the quasiparticle states in the continuum, which also have 
sizable contribution to the pair correlation in the case of nuclei with a shallow Fermi energy 
close to zero. Note here that the quasiparticle states whose excitation energy exceed the 
separation energy form the continuum spectrum because they couple to scattering waves 



10|. The Hartree-Fock single-particle orbits emerge as resonant quasiparticle states with 



finite width 



10|. To describe this situation, one needs to guarantee the asymptotic form of 



the quasiparticle wave function as the scattering wave. In this way, we can avoid artificial 
discretization of the continuum spectrum, and can evaluate the width of a resonant quasi- 
particle state. This allows us to investigate how the resonant quasiparticle states contribute 
to the pair correlation. 

Thirdly, it is important to describe selfconsistently the pair potential, which is the key 
quantity describing the pair correlation. To achieve this, however, the continuum quasi- 
particle states including both resonant and non-resonant states are to be summed up in 
evaluating the one-body densities. We adopt the Green's function technique [l^ that pro- 
vides a simple and effective way of summing. Thus, we are able to perform in the present 
work the fully selfconsistent continuum Hartree-Fock-Bogoliubov calculations, i.e., we derive 
selfconsistently both the Hartree-Fock potential (using the Skyrme functional) and the pair 
potential (using a density-dependent contact interaction as the effective pairing force). It is 
noted that the theoretical framework of the present analysis shares many common aspects 



with that in Refs. 



18N21|. but we differ in that we utilize the self-consistent pair potential 



as well as the Hartree-Fock mean-field. Our ap pro ach is rather similar to the continuum 



Skyrme-HFB calculations formulated in Refs. 



7|,ll3|,l27|l. 



We will perform numerical analysis for the = 86 isotones in the Mo-Sn region. The 
Skyrme-HFB theory with the parameter set SLy4 suggests the presence of weakly-bound 
neutron single-particle orbits above the = 82 shell gap in neutron rich nuclei with ^ 84 
and Z ^50. In the = 86 isotones, particularly, the 3pi/2 and 8^3/2 orbits emerge close to 
the zero energy with the Hartree-Fock single-particle energies ranging from e ~ —0.5 MeV to 
unbound resonances around e ~ 0.7 MeV. Hence it provides us with a good testing ground 
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to study the role of weakly bound low-/ orbits in the pair correlation. 

In the numerical analysis, we shall pay special attentions to the pair density (called also 
the pairing tensor or the abnormal density in the literature) and the related quantities. The 
pair density is one of the most relevant quantities to the pair correlation since it shows up 
in the definitions of both the selfconsistent pair potential and the pair correlation energy. 
Using this quantity we will show that the / = 1 weakly bound or unbound orbits keep finite 
and sizable contribution to the pair correlation even when they form very broad quasipar- 
ticle resonances, and when they are located above the potential barrier. The pairing gap 
associated with the unbound quasiparticle states are also found to stay finite. This result 



is different from those in Refs. 



18 



21|. which suggest the decoupling of the / = 0, 1 weakly 
bound and unbound orbits from the pair correlation. We shall discuss the origin of the 
difference. 

Finally, we remark that the present analysis is related very closely to Ref. [l^, where, 
however, the Hartree-Fock potential is replaced by a simple Woods-Saxon potential although 
the deformation effect is taken into account. In the present work, we do not discuss the de- 
formation effect for simplicity, but instead we investigate in detail the roles of weakly bound 
and unbound orbits especially with the low angular momentum by using the fully selfcon- 
sistent continuum Hartree-Fock-Bogoliubov calculations assuming the spherical symmetry. 

In section II, we describe the formulation of the continuum Skyrme-HFB theory using 
the Green's function technique. After presenting the results including the numerical details 
and the related discussions in section III, we draw conclusions in Section IV. 



II. FORMALISM 

A. Hartree-Fock-Bogoliubov equation with Skyrme interaction 

In the Hartree-Fock-Bogoliubov (HFB) theory, the pair correlated nuclear system is de- 
scribed in terms of the independent quasiparticles. The HFB equation for the quasiparticle 
wave function (piira) in the coordinate space is 

I hira, r'a') - A5(r - r')6^^, h{ra, r'a') V , , ^ ^ , , 

dr y \ ^ (pi{r a ) = Ei(pi{ra) 

^ \ h*{ra, r'a') -h*{ra, r'a') + \S{r - r')6^^> J 

(1) 
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7^(r(T, r'a') 



(2) 



where Ei is the quasiparticle energy, and A is the chemical potential or the Fermi energy. The 
Hartree-Fock Hamiltonian h and the pair Hamiltonian h can be obtained by the variation 
of the total energy functional with respect to the particle density matrix p{ra, r'a') and 
pair density matrix p{ra, r'a') respectively. The two density matrices can be combined in a 
generalized density matrix TZ as 

p{ra, r'a') p{ra, r'a') 

p*{ra, r'a') 6{r - r')5aa' - p*{ra, r'a' 

where the particle density matrix p{ra, r'a') and pair density matrix p{ra, r'a') are just the 
"11" and "12" components of TZ respectively. 

The energy density functional of the Skyrme interaction is constructed with the local 
densities, such as the particle density p{r), the kinetic-energy density T(r), and the spin- 
orbit density J{r), etc., defined with the particle density matrix p{ra,r'a' 
adopt a density dependent delta interaction (DDDl) for the p-p channel: 
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29|. We 



v^^,Xr.r') = ^{l-P^)Vo 



1—7] 



Pir] 
Po 



6{r - r') 



(3) 



which presents similar properties as the pairing interaction with finite range 
pair Hamiltonian h is reduced to the local pair potential 



30j. Then the 



A(r) 



V 



Po 



p{r) 



(4) 



where the local pair density p{r) is defined with the pair density matrix p{ra,ra) j^. 

For the spherical symmetry, the generalized density matrix TZ can be expanded on the 
spinor spherical harmonics as 



7^(ra,rV) = l^,„,(f a)7^,,(r, r')l^*„(fV). 

Ijm 

Using the 11 and 12 components of TZij, one can write the radial local densities as 



Airr 



H\]{r, r') 



J{r) 
p{r) 



dr 



l{l + l) 



T^iHr,r), 



(5) 

(6a) 

(6b) 
(6c) 
(6d) 



where -— denotes the derivative operator with respect to r' acting from right to left. 
dr' 

The quasiparticle wave function is represented as 



rcr 



{r)Yijmircr), where 



which obeys the radial HFB equation 
/ d ff d 



\ 



dr 2m* dr 

A(r 



+ Uijir)-X A(r) 

d fi^ d 



\ 



Uij{r) + X 



(7) 



r,E) = E(Pi,{r,E). (8) 



dr 2m* dr 

The explicit expressions of the effective mass m* and the Hartree-Fock potential Uij{r) can 
be found in Ref . 0, Q , and they are constructed by the radial local densities ([6]) and their 
derivatives. 



B. HFB Green's function and densities with correct asymptotic behavior 

In the conventional Skyrme HFB theory, one solves the radial HFB equation ([8]) with 
the box boundary condition (f)ij{r,E) = at the edge of the box r = R {R being the box 
size) to obtain the discretized eigen solutions for the single-quasiparticle energy and the 
corresponding wave functions. Then the generalized density matrix TZ can be constructed 
by a sum over discretized quasiparticle states. Although the box boundary condition is 
appropriate for the deeply bound states, it is not suitable for the weakly bound and the 
continuum states unless a large enough box size is taken. 

Here the Green's function method is used to impose the correct asymptotic behaviors on 
the wave functions especially for the continuum states, and to calculate the densities. 

The HFB Green's function Qo,ij{f',^' ■, E) can be constructed with solutions of the radial 
HFB equation (IHl). Suppose (f)[Y\r,E) and (f)[^^\r,E) {s = 1,2) are independent solutions 
of the HFB equation ([8]) that satisfy the boundary conditions at the origin and at the edge 
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3l| 



of the box, r = R, respectively, then the HFB Green's function is given by 

go,i,ir, r', E)= J2 ^if(E) [^(^ - E)^lf'^{r\ E) + e{r' - r)0!7')(r, i?)0!f ^^(r', E) 

(9) 



;,s' = l,2 
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The coefficients c^' {E) are expressed in terms of the Wronskians 













H 








14' 







u;/j(r2,+l) wij{r2,+2) 



(10) 



with 



''S('-)i:^i:;''('-)-^s;'(r);^»>i:5(r) 



dr 



(11) 

To impose the correct asymptotic behavior on the wave function for the continuum states, 
we adopt the boundary condition as follows, 



')ij^\r^ E) : regular at the origin r = 



y, E) : outgoing wave at r — )■ oo 



(12) 



Explicitly, the solutions (j)^^\r^E) at r > R satisfy 




r,E) 







(13) 



Here k±{E) = a/ 2m(A ± E)/h and their branch cuts are chosen so that lmk± > is satisfied. 

The generalized density matrix can be obtained by the contour integral of the Green's 
function on the complex quasiparticle energy plane, which in the spherical case can be 
written as [lol, 31 1 

1 r go,i,ir,r',E) ^^^^ 



1 

27n 



dE 



Ce 



The contour Ce should be chosen to enclose the negative energy part of the quasiparticle 
spectra, as shown in Figure [H so that all the quasiparticle states inside the contour are 
summed up. Here the discrete quasiparticle states are denoted by crosses above the Fermi 
energy A. Below the Fermi energy, the continuum quasiparticle states are denoted by the 
solid stripe. As a result, the radial local densities can be calculated by the contour integral 
of the radial Green's function. In this way, we realize a fully selfconsistent continuum Skyrme 
HFB calculations. 
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C. Numerical details 



For the Skyrme interaction, we choose the parameter set SLy4 32], and for the pairing 
interaction the DDDI parameters in Eq. (jl]) are adopted as Vq = —458.4 MeV fm~^, rj = 0.71, 



a = 0.59, and po = 0.08 fm ^ 33l-l35| . They reproduce the experimental neutron pairing 



gap along the Sn isotopic chain. We remark also that the parameter Vq is chosen in such a 
way that the DDDI reproduces the scattering length a = —18.5 fm of the bare nuclear 
force in the low density limit p(r) 0, i.e., in the free space outside the nucleus. Because of 
this constraint and the density- dependence, the pairing interaction strength is large around 
the surface, and even increases in the exterior. The truncation of the quasiparticle states is 
up to the angular momentum I = 12 and j = 25/2 and to the maximal quasiparticle energy 
Ecut = 60 MeV. 

For the contour integration of the Green's function, the path C^; is chosen to be a rect- 
angle as shown in Figure [1] with the height 7 = 1 MeV and the length Ecut = 60 MeV, 
which symmetrically encloses the real negative quasiparticle energy axis. For the contour 
integration we adopt an energy step AE = 0.01 MeV on the contour path. We have checked 
that for the choice of these contour path parameters, the precision for p(r) and p(r) thus ob- 
tained are up to 10~^° fm~^ and 10~^ fm~^ respectively. We choose the box size R = 20 fm, 
and the mesh size Ar = 0.2 fm. We have also checked that dependence of the results on the 
box size is very small thanks to the boundary condition (fT3|) with proper asymptotic form. 



III. RESULTS AND DISCUSSION 



In this section, taking the isotonic chain = 86 as an example, we will discuss in detail 
how characters of the weakly bound and unbound quasiparticle states of neutrons vary as 
the neutron Fermi energy approaches zero, and how they contribute to the pair correlation. 



A. HFB ground states and the quasiparticle spectra 

Some properties of the HFB ground state for ^^^Sn, ^^^Cd, ^^^Pd, ^^°Ru, and '^^^Mo are 
listed in the first four rows in Table HI The neutron Fermi energy (the first row) monotonically 
increases from —2.39 MeV in ^'^^Sn to —0.42 MeV in ^^^Mo as the proton number Z decreases. 
This is because the neutron Hartree-Fock potential becomes shallower as Z decreases, and 
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we could not find a self-bound HFB solution (with A < ) in ^^^Zr and lighter isotones. 
Table [T] also shows the total neutron pair correlation energy 



E ■ 



^-Jdr A(r)p(r), (15) 
and the average pairing gap 

_ Jdr A(r) p(r) 

Jdrpir) ■ ^^'^ 

The quantity in the denominator is the total neutron 'pair number' 

N = p{r)dr, (17) 



which represents the amount of the pair condensate. It is noted that the variation of the 
average pairing gap and the total pair correlation energy along the isotones is less than 10% 
from ^'^^Sn to the last bound nucleus ^^*Mo. 

It is useful to investigate properties of individual quasiparticle states which are the ele- 
mentary mode of single-particle motion in the HFB theory and the building blocks of the 
pair density. It is noted that the spectrum of the quasiparticles, i.e., the eigenstates of the 
HFB equation, includes both the discrete (0 < < |A|) and continuum {E > |A|) quasipar- 
ticle states. Accordingly, the pair density can be expressed as a sum of contributions from 
individual quasiparticle states as 

rEcut 

P{r) = Yl /^-'^•(^) + E / Z^^^ ^^' ^) (18) 

nlj,E„i,<\X\ Ij -^1^1 

where the first term in r.h.s. is the sum over the discrete states, and the second term repre- 
sents the integral of the contribution to the pair density p(r) from the continuum quasipar- 
ticle state with quantum number Ij at energy E . If we include the discrete quasiparticle 
states in the definition of pij{r, E), the above equation can be expressed as 

rEcut 

Pir) = y'pij^'^), where pij{r) = / pij{r,E) dE, (19) 
ij 

and pij{r,E) can be calculated as 

PiAr.E) = (^^±Pj llmgif^{r,r,-E-^e). (20) 
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We can also calculate contributions from the state with quantum number Ij at energy E 
to the pair number as 

hijiE) = J AiTr^pij{r,E)dr, (21) 



which satisfies 



cut 



^ = y2 ni,{E)dE. (22) 

We call the quantity fiij (E) the 'pair number density' in the following. We can also calculate 
the 'occupation number density' nij{E) = J Aur"^ pij{r, E)dr^ which is discussed in Refs. 
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27|. In the following, we will investigate the pair number density nij{E) rather than 



the occupation number density nij{E) since we found that the pair number density nij{E) 
represents more clearly the structure of continuum quasiparticle states relevant to the pair 
correlation. 

With the smoothing parameter e in Eq. fl2U]) the 5 function (no width) originating from 
a discrete quasiparticle state is simulated by a Lorentzian function with the full width at 
half maximum (FWHM) of 2e. In the following calculation, we take e = 5 keV to discuss 
the structure of pair number density. Subtracting the smoothing width 2e = 10 keV from 
the FWHM of the peak, we obtain the physical width F of the peak. 

Figure [2] shows the pair number densities hij (E) for neutron quasiparticle states in a 
low-lying energy interval = ~ 4 MeV in the = 86 isotones. A peak structure below 
the threshold energy E = \\\ (the dashed vertical line) is a discrete quasiparticle state 
(simulated by the Lorentzian function), and a peak above the dashed line may be identified 
as a quasiparticle resonance. The width F of the peak as well as the peak energy -Eq.p. are 
tabulated in Tabled In the same table, we also show the corresponding Hartree-Fock single- 
particle energy e, which is the eigen solution of the Hartree-Fock single-particle Hamiltonian 
h obtained with the box boundary condition. An eigen state with e < is discrete (bound) 
single-particle orbit, and e > is discretized continuum single-particle orbit whose energy 
gives an estimate for the Hartree-Fock single-particle resonance energy. 

In ^^^Sn, the quasiparticle states corresponding to weakly-bound Hartree-Fock single- 
particle states, 2/7/2 and 8^3/2, are discrete states with no width (less than 0.1 keV in 
the actual numerical calculation) as they are located below the threshold energy |A|. The 
Hartree-Fock single-particle state 3pi/2 is already in the continuum [e ~ 70 keV), and forms 
a quasiparticle resonance located just above the threshold energy with i^q.p. ~ |A| + 15 keV 
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with a large width T = 122 keV. On the contrary, the peaks for higher angular momentum 
states, e.g., 2/5/2 and 1/19/2, have finite but smaller widths, forming narrow quasiparticle 
resonances. 

As the proton number decreases, the neutron potential becomes shallower, and both the 
Fermi energy and the single-particle energies are raised up. The width grows larger in ^^^Cd 
than in ^^^Sn for all the quasiparticle resonances above the threshold, although the weakly 
bound single-particle states 2/7/2 and 8^3/2 remain the discrete states in the quasiparticle 
spectrum. 

When the Fermi energy is raised up further in ^^^Pd, the 8^3/2 quasiparticle state becomes 
a resonance located just around the threshold energy, with a large width 95 keV. Since the 
corresponding Hartree-Fock state is still bound (although the binding energy is very small, 
e ~ — 10 keV), it is the pair correlation that makes the corresponding quasiparticle state 
unbound with the finite width. 

In ^3°Ru, the single-particle energies are raised further. All the single-particle orbits under 
discussion except 2/7/2 become unbound. Correspondingly the widths of the quasiparticle 
resonances 8^3/2, Spi/2, 2/5/2 and 1/19/2 grow also. 

When we come to the last bound nucleus ^^*Mo in this isotonic chain, all the quasipar- 
ticle states lie in the continuum. It is noticed that both the 8^3/2 and 8pi/2 quasiparticle 
states become very broad resonances whose width is comparable with the resonance energy 
Eq,p. — |A| measured from the threshold. The quasiparticle state 2/7/2 is now a quasiparticle 
resonance embedded in the continuum although the Hartree-Fock single-particle state 2/7/2 
is still a bound orbit. The width comes from the continuum coupling caused by the pair 
correlation, as the 8^3/2 state in ^^^Pd. The continuum coupling effect is much larger for the 
8p3/2 state than 2/7/2, as is seen in the considerably different values of the width. 

B. Pairing effects on the resonance width 

The large widths of the 8pi/2 and 8^3/2 quasiparticle resonances have two origins. One is 
the barrier penetration of the Hartree-Fock plus centrifugal potential, which is low for the 
states with low angular momentum (the barrier height of the pi/2 and ^3/2 states in ^^^Mo is 
0.851 and 0.844 MeV, respectively). This is effective even without the pair correlation. The 
other is the presence of the pair potential, because of which even a bound single-particle 
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orbit can couple with continuum states, then forms a quasiparticle resonance. The latter 
effect may be seen by examining the pair number density nij{E) under a variation of the 
strength of the effective pairing interaction. 

Table ini shows the quasiparticle resonance widths in ^^®Mo obtained with different pairing 
interaction parameter t] = 0.84,0.71 and 0.62 in Eq. It is seen that for the changes 
of the average pairing gap A,j^ = 0.41 ~ 0.68 ~ 1.08 MeV, the variation of the widths 
of the 3pi/2 and 8^3/2 quasiparticle resonances are F = 0.73 ~ 0.80 ~ 0.94 MeV and 
0.34 ~ 0.40 ~ 0.49 MeV, respectively. We deduce that the pairing effect on the widths are 
approximately ~ 100 keV for 3pi/2, and a slightly smaller for 3p3/2- Comparably large pairing 
effect is also seen in the width of 2/5/2 quasiparticle resonance where the corresponding 
Hartree-Fock single-particle state is a broad resonance already without the pair correlation. 
(On the contrary, the pairing effect is not large for the other quasiparticle resonances arising 
from the hole orbits and those with narrow resonances.) The pair correlation increases 
significantly the width of the quasiparticle resonances corresponding to weakly bound orbits 
with low angular momentum or to the Hartree-Fock single-particle resonances close to the 
barrier top. Since the wave functions of these quasiparticle resonances have significant 
amplitude in the barrier region, the influence of the pair potential on the continuum coupling 
can be effective. 

C. Contribution of continuum quasiparticle states to the pair correlation 

Let us now investigate how the quasiparticle resonances, shown in Fig. |2l contribute to 
the neutron pair correlation. 

For this purpose, we examine their contribution to the neutron pair density in the low- 
lying energy interval = ~ 4 MeV. We denote pnij{r) for the partial contribution from 
the low-lying quasiparticle state, and evaluate it by performing the integral in Eq. ( IT9l) with 
Ecut = 4 MeV for the quasiparticle resonances 2/7/2, 3p3/2, Spi/2 etc. The quantity weighted 
with the volume element, A7rr^p'^ij{r), is shown in Figure |3] for the = 86 isotones. 

As the nucleus becomes more and more weakly bound from ^^^Sn to ^^^Mo, a significant 
variation of the pair density p'nijij") is seen for the 3p3/2 and 3pi/2 states: the amplitude of 
Pnij{r) increases dramatically at the positions far outside the surface, r ^ 7 — 15 fm. For 
the 3pi/2 state, the increase at r = 8 and 10 fm are 80% and 200%, respectively while the 
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increase inside, e.g. at r = 2 fm, is only ~ 20%. The other quasiparticle states 2/7/2,2/5/2 
and 1/19/2 exhibit a similar trend of extending outside but to a much weaker extent. 

Evaluating the volume integral of Pniji^), we list in Table [T] the quantity N^i^j = 
J Airr'^ p'^ij[r) dr, which represents a contribution to the pair number from the low-lying 
quasiparticle states. It is clear that the pair number N^i^ of the 8^3/2 and 3pi/2 states in 
the most weakly bound ^^^Mo is twice as big as those in ^^^Sn, as a result of the obvious 
increase of Pniji^) at the positions r ^ 7 — 15 fm. 

We show also in Table |T] a partial contribution 

4air,n^, = ^ / 4^^'^^ Hr)p'nl,ir) (23) 

to the pair correlation energy from the low- lying quasiparticle state. It can be seen that, 
moving from ^^^Sn to ^^^Mo, the contributions of the 3pi/2 and 3p3/2 states do increase up 
to around 70%. It is clear that this increase of the pair correlation energy is due to not 
only the increase in p'nijij) but also the large spatial overlap between the pair density p'nijij) 
and the pair potential A(r). In Figure H] we show the neutron pair potential A(r) and the 
product of the pair potential and the pair density 4-7rr^A(r)p^^^(r), the integrand of the pair 
correlation energy, Eq. fl23|) . Since the radial profile of the pair potential A(r) is not only 
surface-peaked but also extends up to r ~ 10 fm, the two quantities have significant overlap, 
and hence the product Anr"^ A{r)p'^i-{r) exhibits a significant increase at r 7 — 10 fm. This 
brings about large increase of the pair correlation energy. 

We note here that the energies of the Hartree-Fock single-particle orbits correspond- 
ing to the 3^3/2 and 3pi/2 states move upward around and beyond the threshold: from 
e = —0.48 MeV (weakly bound) to ~ 0.36 MeV (unbound with the single-particle energy 
comparable with the barrier height ~ 0.34 MeV) in the 3p3/2 case, and from ^ 0.07 MeV 
(around the threshold) to ^ 0.68 MeV (even above the barrier height ~ 0.35 MeV) in the 
3pi/2 case. In the least bound case (^^*Mo), the 3p3/2 and 3pi/2 quasiparticle resonances 
have the widths F = 0.396 and 0.802 MeV, respectively, which are comparable with the 
resonance energy. Nevertheless, both the pair number N^ij and the pair correlation energy 
-^pair,zj continue to increase as is seen above. This indicates that the weakly bound and 
unbound states can feel the pair potential and contribute to the pair correlation in sizable 
way. 
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D. Effective pairing gap of continuum quasiparticle states 



In order to make quantitative estimate for the influence of the pair potential on the low- 
lying quasiparticle states, we evaluate state dependent effective pairing gap which can be 
deflned by 

f 47rr^A(r)p' ,.(r) dr 

'-^uv.nlj Pa .2-/ '^^paiv,nlil nlj \^^) 



using the pair density Pnij{r) for the speciflc quasiparticle state. We list it in Table [H 

In the nucleus ^^^Sn, the effective pairing gaps of the 8^3/2 and 3pi/2 states are A'^y^^ij — 
0.591 MeV and 0.570 MeV, respectively, which are about 77 ~ 80% of the total-average 
pairing gap A^°^' = 0.736 MeV. The effective pairing gaps still keep flnite values of the same 
order in the last bound nucleus ^^^Mo, where the 3pi/2 and 3p3/2 Hartree-Fock orbits are 
both unbound. For the 3p3/2 state, the effective pairing gap is A'^^^nij = 0.529 MeV. It stays 
at the level of 78% of the average gap AJ;°*' = 0.678 MeV. Even for the 3pi/2 resonance with 
large width, the effective pairing gap 0.501 MeV keeps 74% of the total. The variation of the 
effective pairing gaps of the 3p states from ^^^Sn to ^^^Mo is also small, i.e., they decrease 
only slightly by ~ 10%. 

The facts that the effective pairing gaps of the 3p3/2 and 3pi/2 states are slightly smaller 
than the total average value, and that they decrease as the orbits become less bound and 



become unbound in the continuum, can be ascribed to the decoupling effect 
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-|21|. which 



is expected to originate from the possible small overlap between the single-particle wave 



function and the pair potential. In Ref. [20|, the effective pairing gap in the weakly bound 
p orbit is suggested to be less than 50% of the average, and possibly less than 1/3 for 
an unbound p orbit. Compared with these numbers, the decoupling effect observed here 
20 — 25%) is much smaller. Namely, we can see that these quasiparticle states persist to 
feel the pair potential and contribute to the pair correlation even if they become unbound 
and have large width. 



The difference between the conclusions of our analysis and those of Ref. [20| can be 
explained as follows. In our analysis, the self-consistent pair potential not only peaks around 
the surface (r 5 — 7 fm), but also extends outside (up to 10 fm or more) as shown in 
Figure H] (a), whereas the pair potential in Ref. 18N21| has a Woods-Saxon shape whose 
main part is concentrated inside the nucleus. Meanwhile, the pair density 47ir'^p'^ij{r) also 
peaks around the surface and extends outsides. Let us take, for instance, the 3pi/2 resonance 
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state in ^^^Mo whose Hartree-Fock single-particle energy is around 0.7 MeV. Considering 
the wave function of the quasiparticle state at the peak energy -Eq.p., its upper component 
(fi{r,Eqp ) has large amplitude around the barrier as the state is located above the barrier 
height (0.35 MeV), and it oscillates in the asymptotic region. On the other hand, the lower 
component (p2{r,Eqp ) exhibits an exponentially decaying asymptotics oc exp(— Kr) with 
K = ^y2m{\X\ + Ec^,p,)/h [l^. Since the contribution of this state to the pair density Pnij{r) 
is given by the product of (pi{r,Eq p) and (p2{r, E^ p), the pair density is confined around 
the nucleus (r < 15 fm in the present numerical examples) and the largest amplitude of 
4:nr'^p'^ij{r) shows up around the surface even though the quasiparticle state is located far 
above the threshold and has a large width. Consequently, such states have sizable overlap 
with the pair potential and thus keep finite effective pairing gap even when the nucleus 
becomes more weakly bound. 

Conversely we may argue a condition for occurrence of the strong decoupling in a semi- 
quantitative way. Since the spatial extension of the pair density is characterized by a size 
constant = 1/k = h/ -\/2m(|A| + -Eq.p.), the strong decoupling can be expected when 
Tp ^ -Rsurf., i-e. only when both the Fermi energy |A| and the quasiparticle energy -Eq.p. are 
sufficiently small. Here the energy -Eq.p. of a discrete or resonance quasiparticle state has a 
lower bound Eq,p, ^ A'^ij given by the effective pairing gap. For the p states in ^^^Mo, we 
find -Eq.p. ~ 0.7 MeV, and |A| > 0.4 MeV, therefore, rp ~ 4 fm, which is not much larger 
than the nuclear radius -Rsurf. ~ 5 fm. This explains why the decoupling is weak here. 

IV. CONCLUSIONS 

We have investigated the neutron pair correlation in neutron-rich nuclei with small neu- 
tron separation energies by means of the fully selfconsistent continuum Skyrme HFB theory, 
in which the Green's function method is utilized to describe precisely the asymptotic be- 
havior of scattering waves for the unbound quasiparticle states in the continuum. We have 
clarified how weakly-bound and unbound neutron orbits contribute to the pair correlation 
properties, especially the orbits with the low angular momentum I = 1 which have large 
spatial extensions. We have chosen the even-even = 86 isotones in the Sn-Mo region for 
numerical analysis, and investigated in detail the pairing properties associated with the neu- 
tron 3p3/2 and 3pi/2 orbits, whose Hartree-Fock single-particle energies (resonances) vary in 
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the interval of —0.5 MeV < e < 0.7 MeV, covering both weakly-bound and unbound cases. 

We found the following features from the numerical analysis. When the 3p quasiparticle 
states are embedded in the continuum above the threshold, they immediately become broad 
resonances with large widths. This is because the barrier height of the Hartree-Fock plus 
centrifugal potential is low for the p orbits (~ 0.35 MeV in the present examples), and also 
because the pair potential which remains effective around the barrier region gives rise to 
additional coupling to the scattering wave in the exterior. The numerical results show that 
the width of the quasiparticle resonances of the 3p states are comparable to the excitation 
energy measured from the threshold. In spite of such large width (F ~ 1 MeV), the contribu- 
tion of the broad quasiparticle resonances to the pair correlation remains finite or can even 
increase. We found that the effective pairing gaps of the broad quasiparticle resonances have 
a comparable size to the total average pairing gap, indicating that the continuum quasipar- 
ticle states persist to contribute to the pair correlation. To be more precise, there exist some 
reduction of the effective pairing gap of 20 — 25% from the total average gap. However, this 
reduction of the effective pair gap is much smaller than what is discussed in Ref. 20|. 

Summarizing, even the broad quasiparticle p-wave resonances in the continuum do con- 
tribute to the pair correlation as far as it is located not far from the Fermi energy. This 



is different from the decoupling scenario 18N21|. The reason for the difference is that the 



pair correlation in the present study is described selfconsistently using the effective pairing 
interaction which has enhancement outside the nuclear surface, and in this case the pair 
potential is enhanced largely around the surface and proximate exterior, keeping overlap 
with the low-/ resonant quasiparticle states. 
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TABLE I: Neutron threshold energy |A|, total average pairing gap A*°*', pair number A^t°*- and pair 
correlation energy -6'*°*^ in the = 86 isotones are listed in the first four rows. The following rows list 
properties of the individual low-lying quasiparticle states shown in Fig. [21 the peak energy i?q.p. and 
the width T extracted from the pair number density hij{E), the Hartree-Fock single-particle energy e 
corresponding to the quasiparticle state, the pair number N^ij, the pair correlation energy -Ep^ir nlj 
effective pairing gap A'^^ evaluated for the quasiparticle state Ij within the energy interval E = — A 
MeV. The unit of the energy is MeV, except for the width F shown in keV. 





i36Sn 


i34Cd 




130Ru 


128Mo 


|A| 


2.390 


1.884 


1.383 


0.894 


0.421 


A tot. 


0.736 


0.721 


0.707 


0.694 


0.678 


jytot. 


16.875 


17.083 


17.458 


18.111 


19.245 


ptot. 
pair 


-6.212 


-6.162 


-6.173 


-6.280 


-6.527 





£ 


-2.302 


-1.799 


-1.297 


-0.801 


-0.309 




^q.p. 


0.760 


0.749 


0.740 


0.734 


0.733 




F 


< 0.1 


0.1 


< 0.1 


2.5 


7.3 


2/7/2 


^uv,nlj 


0.755 


0.745 


0.737 


0.731 


0.727 






3.951 


3.954 


3.959 


3.966 


3.973 




F' 

paiv,nlj 


-1.492 


-1.473 


-1.458 


-1.449 


-1.445 




£ 


-0.480 


-0.235 


-0.010 


0.190 


0.362 






1.965 


1.685 


1.365 


1.081 


0.803 




F 


< 0.1 


< 0.1 


95.0 


238.8 


396.4 


3p3/2 


^uv,nlj 


0.591 


0.579 


0.568 


0.554 


0.529 






0.528 


0.591 


0.679 


0.813 


1.040 




F' 

paiv,nlj 


-0.156 


-0.171 


-0.193 


-0.225 


-0.275 




£ 


0.066 


0.257 


0.423 


0.563 


0.678 




Eq.p. 


2.405 


2.089 


1.760 


1.422 


1.071 




F 


122.3 


264.7 


427.4 


612.2 


802.0 


3pi/2 


^uv,nlj 


0.570 


0.553 


0.540 


0.522 


0.501 






0.179 


0.199 


0.226 


0.268 


0.339 
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TABLE I: (continued) 







i36Sn 


134Cd 


132p^ 




128Mo 




F' 

pair,n(j 


-0.051 


-0.055 


-0.061 


-0.070 


-0.085 




£ 


0.449 


0.862 


1.237 


1.559 


1.816 






2.919 


2.820 


2.694 


2.540 


2.364 




r 


21.7 


50.6 


110.1 


223.8 


431.9 


2/5/2 




0.736 


0.721 


0.708 


0.696 


0.685 






0.690 


0.682 


0.684 


0.701 


0.736 




F' 

pair ,nlj 


-0.254 


-0.246 


-0.242 


-0.244 


-0.252 




£ 


0.566 


1.215 


1.869 


2.522 


3.163 






3.051 


3.186 


3.334 


3.492 


3.661 




r 


0.3 


0.8 


2.2 


7.0 


24.2 


1/19/2 


^uv,nlj 


0.766 


0.756 


0.751 


0.748 


0.747 






1.214 


1.145 


1.079 


1.016 


0.948 




F' 

paiv,nlj 


-0.465 


-0.433 


-0.405 


-0.380 


-0.354 



TABLE II: Dependence of the ground state pair correlation and the quasiparticle properties in Mo on 
the pairing interaction strength. To control the pairing interaction strength, we vary the parameter 77 of 
DDDI in Eq. ([3]) as = 0.84, 0.71, 0.62, 0.56. We list here the threshold energy |A|, the total average 
pairing gap /S.uv, and the total pair correlation energy -Bpair, the peak energy -Bq.p. and the width T of the 
lowest two quasiparticle resonances for Pz/2^Pi/2^ /5/2) aiid /7/2- The unit of energy is MeV, except for the 
width r shown in keV. 



r] 0.84 0.71 0.62 0.56 

|A| 0.363 0.421 0.572 0.778 

AJ;°*- 0.414 0.678 1.077 1.541 

E^°^^ -2.860 -6.527 -13.902 -24.800 

e -0.330 -0.309 -0.292 -0.283 

2/7/2 ^q.p. 0.428 0.733 1.146 1.591 
r 4.0 7.3 13.9 24.6 
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V 


0.84 


0.71 


0.62 


0.56 




e 


-25.566 


—25.571 


—25.562 


—25.534 


1/7/2 


771 

-fcq.p. 


ZD. 204 


25.155 


25.013 






r 


0.4 


0.4 


0.1 


1.6 




e 


0.355 


0.362 


0.369 


0.373 


3^3/2 


i^q.p. 


0.670 


0.803 


1.049 


1.346 




r 


337.7 


396.4 


492.4 


612.6 




e 


—20.175 


—20.151 


—20.114 


-20.078 


2P3/2 


771 

-fcq.p. 


iy.oi6 


iy. /4z 


19.5 /9 


19.3/9 




r 


0.9 


3.1 


9.1 


19.9 






0.673 


0.678 


0.682 


0.685 


3pi/2 


^q.p. 


0.962 


1.070 


1.290 


1.569 




r 


727.4 


802.0 


936.2 


1127.7 




e 


— 18.553 


— 18.528 


— 18.494 


-18.463 




771 


io. iy4 


io.iiy 


i ^.yb 


L 1 . 100 




r 


2.2 


5.6 


13.5 


25.4 






1.817 


1.816 


1.814 


1.811 


^/5/2 


^q.p. 


2.260 


2.365 


2.5y4 


2.899 




r 


392.9 


431.9 


514.4 


660.0 




e 


-20.970 


-21.007 


-21.040 


-21.048 


1/5/2 




20.609 


20.595 


20.498 


20.338 




r 


0.2 


0.5 


2.5 


7.4 
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FIG. 1: (color online) Contour path to perform the integration of the Green's function on the 
complex quasiparticle energy plane. The path is chosen to be a rectangle with the height 7 and 
the length Ecut- The crosses denote the discrete quasiparticle states. The continuum states are 
denoted by the solid stripe below the Fermi energy A. 
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1 2 
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FIG. 2: Pair number densities hij{E) of neutron quasiparticle states with different Ij around the 
threshold energy in the = 86 isotones, obtained with the self-consistent continuum Skyrme HFB 
theory using Green's function method. The thick dashed line denotes the threshold energy |A| for 
the continuum quasiparticle states. The densit^3of state for /7/2 is divided by a factor of 25. 
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FIG. 3: Neutron pair density 47rr^y0^j^.(r) contributed by the low-lying quasiparticle states shown 
in Fig. [2] for the = 86 isotones, where p'^ijif) = dE pij{r,E). The inserts present the 

same density distribution in a log scale. 
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FIG. 4: (a) Neutron pair potential A(r) in the = 86 isotones. (b) Integrand of the pair correlation 
energy, Airr'^ p'^ij{r)A{r), for the 3pi/2 quasiparticle state in the = 86 isotones, where the pair 
density p'^ijir) is the contribution from the 3pi/2 resonant quasiparticle state which is shown in 
Figure [3l 
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